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Using Gaussian integral transform techniques borrowed from functional-integral field theory and 
the replica trick we derive a version of the coherent-potential approximation (CPA) suited for de- 
scribing (i) the diffusive (hopping) motion of classical particles in a random environment and («) 
the vibrational properties of materials with spatially fluctuating elastic coefficients in topologically 
disordered materials. The effective medium in the present version of the CPA is not a lattice but 
a homogeneous and isotropic medium, representing an amorphous material on a mesoscopic scale. 
The transition from a frequency-independent to a frequency-dependent diffusivity (conductivity) is 
shown to correspond to the boson peak in the vibrational model. The anomalous regimes above the 
crossover are governed by a complex, frequency-dependent self energy. The boson peak is shown 
to be stronger for non-Gaussian disorder than for Gaussian disorder. We demonstrate that the 
low-frequency non-analyticity of the off-lattice version of the CPA leads to the correct long-time 
tails of the velocity autocorrelation function in the hopping problem and to low-frequency Rayleigh 
scattering in the wave problem. Furthermore we show that the present version of the CPA is capable 
to treat the percolative aspects of hopping transport adequately. 

PACS numbers: 05.40.-a 63.50.-x 
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I. INTRODUCTION 



The coherent-potential approximation (CPA) is a very 
successful mean-field theory for treating quenched disor- 
der [l|-[5[ and strong correlations in quantum sys- 
tems. It also proved useful for describingclassical diffu- 
sion [Tol - fl2j and vibrational excitations [ill, [lij in disor- 
dered systems. However, until now, The CPA has never 
been applied to the study of topologically disordered sys- 
tems, as all its versions involve a regular crystalline lat- 
tice with additional quenched disorder. As it is, indeed, 
desirable to describe excitations of quenched-disordered 
systems with no crystalline disorder (amorphous met- 
als, glasses, amorphous semiconductors) we present in 
this paper a field-theoretic derivation of the CPA for the 
two formally equivalent problems of diffusion and wave 
propagation in a topologically disordered environment. 
Our effective medium is not a lattice, but a homoge- 
nous and isotropic continuum. Section 2 is devoted to 
the mathematically equivalent problems of diffusion and 
scalar waves with spatially fluctuating diffusivity/elastic 
modulus. For this model our version of the CPA is de- 
rived and solved. From the numerical solutions it is 
demonstrated that the existence of a dc - ac crossover 
in the diffusion problem corresponds to the presence of 
an excess of vibrational states (boson peak) [151 - 120 ] in the 
wave propagation problem. It is shown that the height of 
the boson peak is unrestricted for certain non-Gaussian 
distributions (inverse-power law distribution, log-normal 
distribution) but restricted for a box-shaped and a Gaus- 
sian distribution. For the non-Gaussian distributions a 
scaling relation between the height of the boson peak 
and its frequency position is shown to hold within the 
CPA. For weak disorder the CPA is shown to reduce to 



the self-consistent Born approximation (SCBA), which 
proved very successful to describe the anomalous vibra- 
tional properties of disordered solids [15j-|20j. In sec- 
tion 3 the CPA for vector displacements in an elastic 
continuum with fluctuating shear modulus is presented 
(heterogeneous-elasticity theory [15l420l |). The vector 
theory is shown to possess the same boson-peak features 
as the scalar theory. 



II. DIFFUSION AND SCALAR VIBRATIONS IN 
A DISORDERED ENVIRONMENT 

A. Models and mathematic correspondence 

Let us start our discussion with considering hopping 
transport of electrons or ions in a disordered semiconduc- 
tor, or - equivalently - amon g th e impurities of a doped 
crystalline semiconductor [2ll |22| . Such a motion of par- 
ticles between specific sites i, j in a disordered material 
can be described by a master equation for the probability 
riiit) for being at site i at time t 



d , \ 



,•(*)) 



(1) 



where the transition (hopping) probabilities per unit time 
Wij can depend on the distance ry = |r< — r 3 -| and/or on 
an energy barrier between the sites i and j (repre- 
senting a symmetrized version of diffusion in a landscape 
of states with disordered local energies with phonon- 
assistcd transitions between them [2l|, |22[). The distance 
dependence is supposed to fall off exponentially, so that 
we do not have a long-range model. 
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It has been shown recently [23[ that by a coarse-graining 
procedure such a transport equation can be transformed 
to a diffusion equation on a mesoscopic scale 



9 , ^ 



Vf(r)Vn(r,t). 



(2) 



with a spatially fluctuating diffusion coefficient D(r), 
which is supposed to be a random variable in the 
three-dimensional space with a suitable distribution 
density P[D(r)]. Performing the Laplace transform 
n(r, s) = Jq 00 dte~ st n(r, t) with s = iu + ewe obtain from 
([2]) the diffusion equation in frequency space (disregard- 
ing the t = term) 



sn(r, s) = VD(r)Vn(r, s) 



(3) 



On a macroscopic scale the (quenched) disorder is known 
to lead to a diffusion equation with a space-independent 
but frequency-dependent complex diffusivity D(s) 



sn(r, s) = D(s)V 2 n(r, s). 



(4) 



-D(s) is the Laplace- Transform of the velocity autocorre- 
lation function Z(t) of the moving particle. If the par- 
ticle carries a charge q, this quantity is related to the 
complex, frequency-dependent conductivity a(s) by the 
Ncrnst-Einstein relation 



a(s) 



(5) 



where 



dn/dfi is the derivative of the number of 



carriers with respect to the chemical potential. In de- 
generate quantum systems this quantity is equal to the 
density of electronic states at the Fermi level, in classical 
systems = n/ksT, where T is the temperature. 
Let us now consider a topologically disordered mass- 
spring system in which the masses (which we suppose to 
be equal to unity) are connected by distance-dependent 
force constants, which we call K^. The corresponding 
equation of motion for the scalar displacements of the 
masses at point i is 



dt 



—u i {t) = -Y,K ij (u i {t)-u j (t)) (6) 



The same coarse-graining procedure, which leads from 
(JlJ to {[2]) produces the following stochastic wave equation 



d 2 

^u(r,t) = Vif(r)V«(r,t). 



or in frequency space 
s 2 u(r, s) = §• 



(r,s) = VK (r)Vit(r, s). 



(7) 



(8) 



where now K(r) = v(r) 2 has the meaning of a space 
dependent modulus, which is equal to the square of the 
wave velocity v. K(r) is again a random variable in the 
three-dimensional space and can be identified with D(r) 
in the diffusion problem. 



On a macroscopic scale one deals with an equation of mo- 
tion in frequency space for scalar wave amplitudes u(r,z) 
with a frequency-dependent, complex sound velocity v(s) 



su(r,s) = if(s)V 2 u(r, s). 



(9) 



Here s = s 2 = —uj 2 + ie and K(s) — v 2 (s) is a complex, 
frequency-dependent elastic modulus, which is equal to 
the square of a complex wave velocity. As in optics the 
imaginary part of v(s), v"(u) is related to the disorder- 
induced mean-free path of the waves, ({uj) and to the 
sound- attenuation coefficient T(ui) by [231 ] 



f 



2ujv"(uj) 



1 



£(u) \v(s)\ 2 2\v{5)\ 



(10) 



In a quenched-disordered system, i.e. a medium with 
either spatially fluctuating density or elastic modulus the 
sound attenuation exhibits Rayleigh scattering [23|, [24| 



r(w) oc uj 4 



for 







(11) 



This means that both v(s) and K(s) have a contribution, 
which varies as s 3 / 2 for small frequencies. If we now 
mathematically identify the diffusion coefficient D(s) 
with K(s) we conclude that D(s) has a low-frequency 
non-analytic s 3 / 2 contribution 



AD(s) = D{s) - D(0) oc s 



3/2 



(12) 



which becomes AD(s) oc s d / 2 in d dimensions. This 
non- analytic asymptotics has been proven to hold for 
any quenched-disordered system governed by equations 
of motion of the form JT]), ©, ©, and ©. As D{s) 
is the Laplace transform of the velocity autocorrela- 
tion function Z(t) one obtains the non- analytic long- 
time asymptotics ("long-time tail") Z(t) oc -t( d + 2 )/ 2 . 
This long-time-tail pr operty , which has been known al- 
ready for some time [2514271 ] . is obviously equivalent to 
Ra y lei g h scattering via the correspondence D(s) <H> K(s) 
[23, [28j. The analogy, in fact, goes further: In a 
quenched-disordered system (which is the subject-matter 
of the present work) it is known that a cross-over hap- 
pens between a frequency- independent diffusivity and a 
strong frequency dependence, which can be parametrized 
as D(s) oc s x with x w 0.8 [IMlj]. This disorder- 
induced diffusivity transforms under the corresondence 
D(s) O K(s) to disorder-induced anomalous frequency 
dependence of the elastic modulus ; the onset of which 
corresponds to the boson peak [H, HH . This correspon- 
dence will be discussed in more detail below. 



B. Derivation of the CPA 

We now consider the mathematically equivalent problems 
Eqs. ^ and ©, identifying the quantities D(s), K(s) 
and s, s. We define [67| 

(s - VL>(r)V)<5(r - r') =< r'\A[D(r)]\r > (13) 
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The Green's function corresponding to Eq. ([2]) is given 
by the inverse matrix element of A: 



G(r,r')=<r'|^ 1 |r> . 



(14) 



Applying standard methods in replica field theory |35[ 
we represent the Greens function as a functional integral 
over mutually complex conjugate fields u(r) a and u(r) a 
present in a = 1, . . . , n replicas of the system as follows: 

G(r,r') = 11/ ■D[u a (r) ) u a (v)}u 1 (r)u 1 (r') 



a=l ' 

-E < u a \A\u a > 

x e a 



(15) 
(16) 



with the generating functional 



« r -J2<u a \A\u a > 

Z[J(r,r')] = [] / V[u a (v)]V[u a (v)}e « 

a=l 

-E < w (q) |J (q) I« (q) > 

e « (17) 

and the source- field J a . By an integration by part the 
matrix elements of the inverse Green operator A can be 
written as 

< u a \A[D]\u a >= J d 3 r(s\u a (r)\ 2 + D(r)|Vw Q (r) 

n (18) 

We now apply the Fadeev-Popov procedure [36[, which 
consists of the replacement of the fluctuating diffusivity 
-D(r) by a complex auxiliary field (r, s) with the help 
of a delta functional, which, in turn, is represented by 
another auxiliary field A^ (r, s) 

Z[J] = Jv[u,u] J V[Q]e- <u \ A[Q] - Jlu> 5[D -Q] 
= f V[u,u]V[Q,A]e-<^ A ^- J ^ > e <A \ D - Q> 



V[Q,A] e 



-Tt{ln[A[Q]-J] } P <A\D-Q> 



(19) 



where we have suppressed the replica indices for brevity. 
The third equality in Eq. (Q~9|) follows from integrating 
out the displacement fields u a and u a . In order to pro- 
ceed further we devise another coarse-graining procedure. 
We tile the total space into N c cells of (approximate) vol- 
ume V c = V/N c , where V — L 3 [68[ is the total volume. 
This could just be done by means of a cubic grid. How- 
ever, in order to avoid any relation to a crystalline lat- 
tice we think, instead, of a a Voronoi tessellation around 
midpoints of a closed-packed hard-sphere structure. This 
gives V c — L z c = (Tr/6)r] c dl, where r] c ~ 0.56 is the close- 
packed packing fraction, resulting in L c w 0.66<i c . Within 



a cell with label i we replace the diffusivity by their av- 
erage in each cell and assume that a diffusion equation 



d_ 

at 



n(r,t) = VAVn(r,i), 



(20) 



holds within a cell with label i. We now assume that the 
random numbers Di arc independent of each other, i.e. 
the joint distribution density is assumed to factorize as 
P(Di . ..D Nc ) = T] p(A). 

i 

Our assumption of independent fluctuations of the quan- 
tities Di implies that the size of the cells L c must be 
larger or at least equal to the correlation length £ of the 
diffusivity fluctuations, which is defined by 



^ = W)j rf HD(r + r )£(r )). 



(21) 



Correspondingly we confine the k summations in the sub- 
sequent analysis to remain below a cutoff |k| < = v/£, 
where v is an adjustable number of the order of 1. 
Within our model D(r) is now a piecewise constant func- 
tion in real space and the same should hold for the auxil- 
iary fields Q and A, which are now labeled as Q^, A\ a \ 
Using this the scalar product, which appears in the ex- 
ponential in Eq. (|19l) can be written as: 



< A\D-Q>=^Y,E A * Q V) (D^ - Q\ a) ) (22) 



We now start to evaluate the configurational average. 
Due to the Fadeev-Popov transformation the only term 
to be averaged over is the term e < A P-Q>. 
Assuming that all the N c coarse-graining cubes behave 
the same on average and using that the individual cubes 
are not correlated, we can write 



<a\d-q>^ = rj Q rL\ e ~ 1 1 ~ * 

= p E a ^m((exp[-^AW(^-QW)] >.) 



(23) 



Note that the two occurring volume ratios do not cancel 
each other due to the average inside the logarithm. Using 
the generating functional (|T^1) can be written as 



Z[j] = / V[Q,A]e 



-S eS [Q,A,J] 



(24) 



where we have now replaced the source field J a (r, r') by 
translational-invariant source field J(r — r') which is not 
supposed to depend on the replica index a. The effective 
action takes the form 



S eS [Q,A,J] 



Tr{ In (A[Q]-J)} 



(25) 



Vr 
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Since the factor y- in the effective action (|25|) is much 
larger than unity a saddle point approximation can be 
employed to evaluate the integral in (f24f . In general this 
factor will scale as 



CO 



(26) 



v_ = (L\ a d 

Vc \U 

Accordingly the CPA becomes exact for d — > oo Q ■ 
We now assume that the fields Q and A are replica inde 
pendent. 



This corresponds to the introduction of an effective ho- 
mogeneous medium. In this medium (I3I[) becomes: 



Vc A = Vcd_ 
V V dQ 



V k s + Qk 



(32) 



In the second step the effective-medium operator 



A eS (k, k) = (s + Qk 2 )5 k ~ k (33) 



S cS [Q, A, 0] = n S' eS ({Qi}, {Ai}) (27a) 



SUiQihi^i}) = tr{m[!(Q)]} 



(27b) 



g-J^A^fl.-Qi) 



was defined. From this representation one can see that 
in the CPA the following holds: 



(G)(k,k,s) 



s + Qk 2 kk 



Dk 2 



( 34 ) 



where "tr" now means a trace whithout the replica in- 
dices. 

The saddle point is determined by the equations 



dS' eB 



dQi 
dS' eS 



OA,, 



Qi=Qi,s 



Ai=Ai 



0, Vt 
0, Vi 



(28a) 
(28b) 



The derivative with respect to A is easily performed and 
yields: 



= 



&i Qi,s \ 



exp[^A iiS (A - Qi, s )] 



(29) 



Since ^ <C 1 the exponential in the denominator can be 
expanded to first order |69l |: 



A - Qt 



1 + ^(A - Qi, s )Ai 



(30) 



The second saddle point equation gives 



dtr{ln[A(Q)}} 



dQ 1 



$A hS (e-^ D > 



e -^-A i:t {Di-Qi 



Vc 

V 



A; 



(31) 



The left-hand side can be evaluated under the assump- 
tion that the saddle point field Qs is constant in space 

Qi,s = Qi Vi 



under the assumption that the averaged system exhibits 
translational invariance. This equation expresses the av- 
eraged Green's function in terms of the Green's function 
of a homogeneous medium, where the spatially fluctuat- 
ing diffusivity is replaced by the self energy Q, which, 
however, is now frequency dependent: The space de- 
pendence due to the disorder has been transformed to 
a disorder-induced frequency dependence. 
From (|32t it follows that if Q is homogeneous in space, 
the same holds for A. Defining a new field A = 3V c /j/VA 
with v = v 3 /2ir 2 and performing the summation in (|32[) 
with a cutoff |fc| < k^ the CPA equations become: 



= 



A - QO) 



1 + jf[A-Q(*)]A(«), 

~ 3 f k z k 2 

A(s) = ± / dkk 2 

fcf Jo s + k 2 Q(s) 



(35a) 



(35b) 



1 



■[l-sG( S )] 



Q(s)' 

with the local Green's function 



3 f 5 I 

G(s) = -3 / dkk 2 

k 3 J s + k 2 



Q(s) 



(35c) 



We call A(s) the susceptibility function, because it is pro- 
portional to the local dynamic susceptibility of the dif- 
fusing particle. 

The CPA equation (|35a|) can be cast into the following 
equivalent forms 



Q(s) 



l + f[A-Q(s)]A( s ), 

A ' 

l + |[A-Q(s)]A( s ), 



(35d) 
(35e) 
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It is worthwhile to note that the k integral in Eq. (|35b|) 
for the susceptibility function can be carried out ana- 
lytically. The diffusion pole in the denominator of the 
integrand produces a A(s) = A(0) + const x s 2 / 3 low- 
frequency asymptotics, which is inherited by the function 
Q{s). So, independent of the type of disorder, we obtain 
a correct long-time behaviour for Z(t) and Rayleigh scat- 
tering for the wave problem. Because for the vibrational 
problem in the s — > limit the disorder scattering is sup- 
pressed by the Rayleigh frequency dependence, the CPA 
expression for the imaginary part of Q(z) can be shown 
to reduce to the Born approximation in agreement with 
the previous derivations [23|, [24j • 

In the dc limit s — Owe have A(0) = 1/Q(s), so we obtain 
from Eq. (|35cl) 



Q(o) 



D, 



V D, 
3 3Q(0) 



(36a) 



(From now on we suppress the index i, which indicates 
the average over p(Di).) In the case Q(0) ^ (which 
is not trivial, see the paragraph on percolation) one can 
divide by 3Q(0)/V to obtain 

(36b) 



C. Relation with previous effective-medium 
theories 



1. Lattice CPA 

The standard 2-site coherent-potential approximation for 
the hopping problem of Eq. (TTJ on a lattice is 



Wij - r( s ) 



i 



+ (W ij -T(s)^(l-sG ii (s) i 



= 



(37) 



Here T(s) is the effective frequency-dependent hopping 
rate and Gu(s) is the local Green's function of the effec- 
tive medium, which whithin this theory is a simple-cubic 
lattice with coordination number Z = 2d = 6. The lat- 
tice Green's function has the form 



Gu(s) = J2 



1 



keBZ 



fl + r(*)/(k) 



(38) 



where the sum goes over the first Brillouin zone (BZ), 
and 

/(k) = 6 — 2[cos(fc :E a) + cos(k y a) + cos(fc z a)] (39) 

with a being the lattice constant. Defining the local sus- 
ceptibility function 



(40) 



the CPA equation (|37l) takes the form 



- r( s ) 



= o 



(41) 



i + (Wy -r( s )iA«( s ), 

Now we take the continuum limit by replacing the BZ 
k summation by ^ k — > ^ dk and /(k) by its low- 

wavenumber limit k 2 a 2 . If we now define the local diffu- 
sivities by Di = W^a 2 , the effective- medium diffusivity 
by Q(s) = r(s)a 2 and the continuum susceptibility func- 
tion by A(s) = Au(s)/a 2 we arrive at the continuum- 
off-lattice CPA result (|35a|) . provided we take v = 1. In 
order to be consistent with the continuum limit of the lat- 
tice CPA one may take always this value. On the other 
hand, one can also use this value such that the CPA per- 
colation threshold p^ PA =V/3 (see below) agrees to the 
continuum percolation threshold p c of a certain topology 

MM. 



2. Self-consistent Born approximation, SCBA 

If one takes Gaussian disorder for the local diffusivity 
one can perform the disorder average over the generat- 
ing functional (fl"7|) exactly, which leads to an interact- 
ing effective field theory with the variance of D(r) as 
coupling constant. Taking apart this interaction by a 
Hubbard-Stratonovich approximation and then perform- 
ing a saddle-point approximation (assuming a small rel- 
ative variance of D(r)) one arrives at the self- consistent 
Born approximation for the scalar problem |4ll443j . We 
can, however recover the SCBA from the CPA in the 
following way. Defining Do to be the average of the fluc- 
tuating diffusivities and defining the quantities Q(s) = 
Dq — S(s), Di = Do — Aj, we obtain from (|35al) the two 
(equivalent) CPA equations 



= 



Aj - E(s) 



£(*) 



l-f(A i -S( S ))A(s) 



1-^(A,-S( S ))A( S ) 



(42a) 



(42b) 



We now expand the interior of the average in (|42b|) with 
respect to Aj — E(s) to lowest nonvanishing order (re- 
specting (Aj) = 0) we obtain 



£(*) 



(A 2 )|A( S ) 



(43) 



which is the SCBA for the scalar problem. As indicated 
already above, the SCBA can also be obtained from the 
saddle-point equation (1291) by putting the exponential not 
into the denominator but into the numerator and then 
expand with respect to the small number V c /V to first 
order. Because then only the first two cumulants of the 
distribution of the Di enter, this corresponds to assuming 
Gaussian disorder. 

So we recover the SCBA from the CPA in the Gaussian 
and weak-disorder limit. 
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3. Network effective-medium approximation, EMA 

The CPA-like effective-medium treatment of the 
impedances of a heterogeneous medium or network date 
back to Bruggeman [44| and Landauer 45 1. For a disor- 
dered if-fold coordinated network of fluctuating conduc- 
tances <?i the expression for the effective-medium conduc- 
tance g m is [HI 







9i + (f - l)ffr, 



(44) 



It has been generalized for the ac problem [3ll . l47j . setting 
Z/2 = d 







9m (s) 



,g% + (d- l)flm(s) + ds 
which can be re-arranged as 



= 



9 - 9m{s) 



l+[<?-5m(s)]3i^ 



(45) 



(46) 



which has - for d = 3 - the same form as the CPA equation 
(135a]) with 



A( S ) 



1 



EMA 



Q(s) 



(47) 



This EMA has the same analytical structure as EMA ver- 
sions derived earlier in the literature [HI, HO, |48| . While 
these ac effective-medium theories describe rather nicely 
measured hopping conductivity data they violate the 
non-analyticity requirement (TT2"1) . As stated above, in 
CPA the function A(s) has a contribution, which varies 
as s 3 / 2 , whereas A(s) E ma does not, as can be clearly seen 
from Eq. (gTJ). 

D. CPA results for the dc diffusivity 

1. Percolation 

In order to treat the percolation problem, which can be 
considered as the continuum version of the Lorentz prob- 
lem [2(| [2?i EH, we assume a distribution of local 
diffusivities of the form 

P(Di) = pS(Di - D ) + (1 - p)5(Di) (48) 

where p is the volume fraction in which the diffusivity is 
non-zero. Inserting this into Eq. (|36a[) one obtains 



Q(0)g =Q(0) 



1 



i+(3 _i)QM 



(49) 



which has, like Eq. (I36a[) always the trivial solution 
Q(0) = 0. For the case Q(0) ^ we obtain from Eq. 

ma 



o(o) = gA>(p-Pc)/(i-p c ) 



(50) 



with p c = ^. For p 
Q(0) = takes over. 



< p c the trivial solution of (|49l) 



2. Activated diffusion 

In this class of models the diffusion of a particle is con- 
sidered to take place by jumps over barriers of height Ei 
with a certain distribution P(e) ("random barrier model" 
13111'). So we write 



Di = D e E ' /kBT 



(51) 



and parametrize the dc diffusivity as (S — 1)Q(0) 
D e E ^ kBT . Then Eq. (j36b| takes the form 



V - = J dEP(E) 



1 



p (E-E a )/k B T + 1 



(52) 



In the low-temperature limit the Fermi function becomes 
a step function and we obtain 



: Pc 



dEP(E) 



(53) 



which means that the parameter E a becomes tempera- 
ture independent. From this follows that in the random- 
barrier problem the dc diffusivity is always of Arrhc- 
nius form, as observed frequently in fast-ion conducting 
glasses [511 . 

It is worthwhile to point out that (153[) corresponds to 
the so-called percolation construction for obtaining the dc 
conductivity of a disordered hopping-conduction network 
(22I l3~fl ] . It has been nicely demonstrated recently that 
the low-temperature physics of the random-barrier model 
is essentially percolation physics [52[. 



3. Variable-range hopping 

A rather widely investigated type of carrier diffusion 
in disordered materials is that of electrons perform- 
ing phonon-assisted tunneling transitions between local- 
ized states ("hopping transport" [HI)- Understanding 
the mechanims of electronic hopping transport has been 
shown recently to be of extreme importance for devis- 
ing organic light-emitting diodes (OLEDs) Here we 
show that by the CPA one recovers the classical results 
of Mott [HI and Efros, Shklovskii [22j for variable-range 
hopping. 

The local diffusivity depends on an activation barrier E 
and a characteristic hopping distance r with distributions 
P{E) and P{r) 



Di = D e 



-ar-PE 



(54) 



where e ar is the tunneling factor and [3 = l/fc^T. De- 
pending on the density of localized electronic states near 
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FIG. 1: Comparison of ionic ac conductivity data of sodium 
trisilicate glass [55|] with the CPA prediction for activated 
hopping with a constant-barrier distribution. We use units 
in which Do = k$ = 1. The arrows indicate the boson-peak 
positions ojbp = oJbp of Fig. [3] 



the Fermi energy P(E) is either considered to be con- 
stant (Mott hopping) or proportional to E 2 (Coulomb- 
gap, Efros-Shklovski hopping). The distribution of sites 



P(r) = -^pr 2 0(R 

z 



R :i 



r 2 9(R - r) 



(55) 



where Z = ^AnpR 3 is the number of adjacent sites within 
a given radius R. If we parametrize the dc diffusivity as 
- 1)Q(0) = -D e" ? , we obtain from Eq. (l36b| 



Pc = J dEP(E) J drP{r) 



1 



e ar+/3E-i + 1 



(56) 



In the low-temperature limit the Fermi function be- 
comes again a step function, and, by means of integra- 
tions by part one obtains the famous results lnQ(0) oc 
-(T /T) 1 /4 (Mott hopping) lnQ(0) oc -(Tq/T) 1 / 2 
(Efros-Sklovskii hopping) . Again, these results are equiv- 
alent to the percolation construction (22[. 



E. CPA results for the ac diffusivity 



We consider activated transport of the form (f5Tj) with a 
constant barrier distribution 



P(E) = — < E < E* 
K ' E* ~ 



(57) 



which is equivalent to an inverse- Power distribution for D 
1 1 



P(D) = 



In fi/a D 



[i < D < a 



(58) 



with a = D and /! = DqC~^ e . 

In Fig. [T] we show ac conductivity data collected from 
the literature over a very wide range of frequencies by 
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FIG. 2: Loss function [Q'(cj) - Q(0)]/w oc [a(u) - ff(0)]/w 
calculated from the CPA curves in Fig. [T] corresponding to 
the temperatures 456 K, 756 K and 1273 K. Below the dc -ac 
crossover the Rayleigh-type non-analyticity is visible. 



Wong and Angell [55j together with the CPA prediction 
for the constant-barrier model (|57[) . The only input is 
the measured dc activation energy of E a = 75 KJ/mole 
frol | However, it should be made clear that the old EM A 
theories [H, H(| IH, H3, Hl| are also able to produce such 
a fit. The difference to the CPA can be seen from the 
"loss function" [HlHi [a(u) - cr(0)]/w oc [Q'(w) - 
which is shown in Fig. [2] At low frequencies this function 
behaves as w 1 / 2 due to the Rayleigh-type non-analyticity. 
The old EMA theories do not exhibit this non-analyticity. 
In the figure this behaviour is demonstrated. There is 
experimental [59| and simulational [6(| evidence for this 
low-frequency non-analyticity of hopping transport. 



F. Model calculations for the scalar phonon 
problem 

We shall now exploit the mathematical correspondence 
ico <H- -lo 2 , Q(s) = D(s) o Q(s) = K(s) and, corre- 
spondingly, discuss the vibrational anomalies induced by 
the quenched disorder, as given in CPA. Similar discus- 
sions have already published in the literature [32| , where 
it was pointed out that the boson peak (BP) in the vi- 
brational problem corresponds to the dc - ac crossover in 
the diffusion problem. 

As shown in the last subsection, the activated-diffusion 
model with a constant barrier distribution corresponds to 
an inverse-power law distribution for the local diffusivi- 
ties. For the corresponding local moduli K we re-write 
this distribution 



P(K) = 



1 1 

In [i I a K 



fi < K < a 



(59) 



Other possible distributions, namely a uniform distribu- 
tion, a Gaussian distribution, truncated at K = and a 
log-normal distribution are detailed in table HI The den- 



sity of states (DOS) of the scalar phonons can be calcu- 
lated as 



2u> 



ffM - — Im{G(s)} 



(60) 



where G(s) is given by Eq. (|35cp . and we identify the 
correlation cutoff k^ with the Debye cutoff hp. 
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Log-Normal 
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FIG. 3: Reduced density of states g(£u) / go(fi>) for inverse- 
power law disorder with the same disorder parameters as in 
Figs. [l]and[2] As the disorder increases so does the height of 
the boson peak, it is also shifted to lower frequencies. 



TABLE I: Probability distributions used to model the disor- 
der. 

As mentioned above, the CPA makes it possible to de- 
scribe highly disordered systems, i.e. systems in which 
the variance of the spatial fluctuations of the quantity 
of interest exceeds the square of its average. In order to 
quantify the strength of the disorder, we define a disor- 
der parameter as the ratio between the variance and the 
squared mean of the disorder distribution 7 = (K 2 )/{K) 2 
For the truncated Gaussian and the uniform distribution 
7 has an upper bound. For the uniform distribution a 
maximum disorder strength of 7 = i can be reached, for 
the truncated Gaussian this limit is 7 = ^ — 1. Thus 
two of the four distributions can only model medium to 
weak disorder. On the other hand, the inverse-power 
and log-normal distributions have no upper bound of the 
disorder parameter. In particular, for the inverse-power 
distribution with p = ae~ fiE * the relation 7 = f3E* /2 
holds. 

In Fig. [3]we show the so-called reduced DOS g(uj)/g£>(Q), 
calculated for the inverse-power-law distribution. We use 
units, in which k^ = fcrj = 1 and Kq = 1. <?d(w) = 
3w 2 /wfj is the Debye DOS, and u D = ^Q(0) is the De- 
bye frequency. 

The disorder parameters 7 = (3E* /2 have been chosen to 
agree to those in the conductivity calculations of Fig. [T] 
The boson peaks shown in Fig. [3] increase with increas- 
ing disorder, while its position decreases. The positions 
wbp of the boson peaks are indicated in Figs. [T]and[2]via 
the correspondence wbp <^bp- ^ is clearly seen that 
the boson peak marks the onset of the disorder-induced 
frequency dependence of the diffusivity (= conductiv- 
ity). Via the correspondence D(uj) O K(Cj) = v(Q) 2 this 
means that the boson peak marks the frequency depen- 
dence of the sound velocity in the disordered vibration 



model. This is in agreement with earlier conclusions from 
effective-medium calculations using the EMA 32, 33] and 
the SCBA 

In Fig. U] the boson-peak height is plotted against its 
frequency position for all four distributions considered. 
We include also the prediction of the self-consistent Born 
approximation (SCBA). Fig. [5] shows a close up of the 
small-disorder region. 

The down-shift and reinforcement of the boson peak with 
increasing disorder has been shown in the SCBA-based 
earlier treatments of vibrational anomalies 20] to re- 
sult from the disorder-induced level repulsion, which is 
present in the anomalous frequency regime above the 
boson peak. This level repulsion, which is typical for 
random-matrix spectra, results from the absence of sym- 
metries on the microscopic scale. 

From our CPA calculations, displayed in Fig. |4] it fol- 
lows that the height of the BP scales with its frequency 
position wbp as f^rj^k ( w bp) c - with c = 1.25. It 
is suggestive that this scaling should be related to the 
power-law frequency dependence of the diffusivity in the 
mathematically equivalent diffusion problem, but we did 
not find a way to prove this. 



III. CPA FOR HETEROGENEOUS-ELASTICITY 
THEORY 

A. Model 

We start with the equations of motion of elasticity the- 
ory in frequency space, formulated in terms of the stress 
tensor aij — p(Jij (p is the mass density) |6lj |: 
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FIG. 4: Relation between boson peak height and position 
for the CPA solutions of different distributions and disorder 
strength for k^ — kr>- 
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FIG. 5: Close up of the low disorder region of Fig. [4] Note 
that in the low-disorder limit ujbp — > 0.5333, independently 
of of the chosen disorder distribution. 



0J 2 Ui 



(r,w) -^djVijfau}) ^^AijUjfau), (61) 



where Uj(r, id) are the Cartesian components of the dis- 
placement field. We consider an elastic medium in which 
the elastic shear modulus may fluctuate in space. If the 
system is assumed to be still isotropic, the stress tensor 
can be represented in two ways 



1 



(62a) 



<7jj = -cry = XSijtT{e} + 2G(r)e y - 

= -cry = KS i:j tT{e} + 2G(r)% (62b) 

Here A = pX = AT+|G is the longitudinal Lame modulus, 
K = pK is the bulk modulus and G — pG the shear 
modulus. e%j is the strain tensor (dj = d/dxj) 



l/2(d i u J + djUi) 



(63) 



-Sijtr{e} 



(64) 



(dj = d/dxj). The spatial fluctuations of the shear mod- 
ulus can be modeled in two ways: In what we call Model 
I. 0, Htjj the longitudinal Lame modulus A is assumed 
to be constant (referring to the representation (|62aj) ), in 
Model II. [17J, [l8| the K modulus (referring to the rep- 
resentation (|62bl) ) is assumed to be constant. As can 
be shown easily, in model I. the macroscopic longitudi- 
nal Lame modulus is frequency independent, in model II. 
the macroscopic bulk modulus is frequency independent 

Ea Hani- 



B. Derivation of the CPA for heterogeneous 
elasticity 

For model I. the matrix Ay takes the explicit form 



At 



-u> Sij — Xdidj 



(65) 



djG(r)di + S tJ d t G(r)d t 



for model II. we have 

Aij = -id 2 5 l3 - Kdidj (66) 

+ ld l G{r)8 J - (d j G{T)d i + 8 ij Y,d t G{v)d l 
6 v t 

The Green matrix G = A^ 1 is represented as 

n . 

G(r,r% = J] V[u^(v),u a m (v)]u\(v)u)(v') 

-E<«?I4<> 

x e « (67) 
Z[J(r,r')] (68) 



with the generating functional 



Z[J(r,r')] = J] / Z>[«?(r)<(r)]e 



a = l 



E < U ? \Jhn\ U m > 



(69) 



and the source-field matrix J a . By an integration by part 
one arrives at the following representation of the action 

J2<<\AlG]\u«> = /'d 3 r^|<(r)| 2 + ^trK(r)} 2 

ij ^ i 

+ G(r)J2\^(r)\ 2 ) (70) 



Applying again the Fadeev-Popov procedure and per- 
forming all the steps we have done before, we arrive at 
an effective action, which looks similar to that of the dif- 
fusion problem (|2~51 



S eS [Q,A,J] = Tr{ In (A[Q] - J)} 



(71) 



n v 



e v 



But now the trace operation goes also over the Cartesian 
indices. 

The homogeneus Matrix ^4 c ff[Q] is both diagonal in the 
Cartesian indices and with respect to the k vectors. The 
diagonal elements arc 1/Gi(k, s) and twice l/GT(k, s), 
which are given by 

l/G L (k,.§) = ~s + k 2 [\ + 2Q(~s)] Model I (72a) 



l/G L (k,S) = S + k 2 [K+ -Q(s)] Model II (72b) 




0.15 
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FIG. 6: Boson peak data for SiC>2 [63J compared with the 
reduced DOS for the inverse-power model with two different 
lower cutoffs fi. 



l/G T (k, 5) = S + fc 2 Q(i) 



(72c) 



The saddle-point equations, followed by the expansion 
of the exponential in the denominator leads to the CPA 
equations 



Gi - Q(S) 



l + H(G i -Q(s))A(§) 



(73a) 



The susceptibility functions for the models a = I, II take 
the form 

A Q (S) = A jf 5 dkk* (q a G L {k, s) + 2G T (fc, s)) (73b) 

with qi = 2 and qn = 4/3. The density of states is 
calculated from the well-known formula 



(/UM = Tm< { +2Gvl«) 



(74) 



where we have introduced the local longitudinal and 
transverse functions (identifying again £ with Ud) 



Gl,t{s) = 4- / ° dkG LiT (k,s) 



(75) 



We can write the susceptibility function A Q (s) as follows 

(76) 



l-SGr(g) , 1-.§G L (S) 



Q(s) ' *"p a + q a Q(s) 

q a Q(s) l-SG L (s) 
Pa + q a Q(s) 1 - sG T (s 



A(5) 1 + 



with pj = X and p// = AT. A(s), which is the transverse 
local susceptibility, is the same mathematical function 



of Q as the susceptibility function of Eq. P5bl) for the 
scalar phonon problem. Because the function inside the 
big brackets is only weakly frequency dependent and the 
density of states is dominated by the transverse Green's 
function all the results derived and presented for the 
scalar phonon problem hold also for the vector phonon 
problem. 

In particular, for non-Gaussian distributions of the shear 
modulus the quantity g(oj)/go(oj) can have boson peaks 
with arbitrary heights, and a scaling as depicted in Fig. @] 
holds [HI]. As an example we show in Fig. [S]the reduced 
density of states, extracted from inelastic X-ray measure- 
ments by Baldi et al. [63[ together with the correspond- 
ing quantity calculated for the inverse-power distribution 
P(G) given by Eq. (53]) 

The important difference between the scalar model and 
the vector theory (heterogeneous elasticity theory) is that 
it describes the physically relevant vector displacements, 
in which dilatational and shear degrees of freedom can 
be distinguished. Model II. represents a theory, in which 
the disorder-induced anomalous frequency dependence is 
dominated by the dilatation-free shear degrees of freedom 
in agreement with recent computer simulations [l7l . [Til . 

Mm- 



IV. CONCLUSIONS 

We have derived a version of the coherent-potential ap- 
proximation for both diffusional and vibrational motion 
in a quenched-disordered environment, which is suit- 
able for topologically disordered materials. The effective 
medium is not a crystalline lattice but a homogeneous 
and isotropic system with frequency-dependent diffusiv- 
ity or elastic constants, resp. . The results can be di- 
rectly applied to experimentally measured spectra. In the 
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weak-disorder limit the CPA has been shown to reduce 
to the self-consistent Born approximation, which is based 
on Gaussian disorder. In the strong-disorder limit, in 
which the local diffusivities or elastic quantities vary ex- 
ponentially, the CPA has been show to correctly desribe 
the percolative aspects of such systems. The disorder- 
induced vibrational anomalies have been shown to be- 
come stronger as the disorder is increased. In particular 
the height of the boson peak has been shown to increase 
indefinitely with the disorder. 

In contrast to earlier effective medium theories for the 
diffusion problem the present theory includes the correct 
low-frequency non-analyticity, which leads to a long-time 
tail of the velocity autocorrelation function of the diffu- 
sion problem and to Rayleigh scattering in the vibra- 
tional problem. 
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From here on all steps are analogous to the previous sec- 
tion. First the replica trick is performed and a coarse- 
grained potential 

V(x) = (80) 

i 

with no correlation between the coarse-graining boxes is 
introduced. V(x) is subsequently replaced in ,4 via the 
Fadeev-Popov precedure by an auxiliary field Q yielding 

(C[J})= J V[Q,N\e~ nS ^^ (81) 
with the effective action 

S* cff [Q,A,J] = -trln(A + J) -Y / ^{^ Vc/VMvi ~ Qi) ) . 

(82) 

As described in the body of the paper, the CPA equations 
determine the saddle point of (l81j) . The saddle point 
equations read: 



Appendix:Electrons in a random potential 

The method developed above for the diffusion and scalar 
phonon problem can, with slight modifications, be ap- 
plied to the problem of an electron gas in a random po- 
tential. This model is governed by the Schrddinger equa- 
tion 



/d»**t(*)(g 



V(x) 
E = 



E) ^(x)\ip >= A\(p >= 



(77) 



E+ie 
N 



The corresponding Green's function can then be ex- 
pressed with a coherent-state path integral over the 
Grassmann fields 9 and 9: 



G(x,x',E) 



1 



6C[J] 



C[0] SJ{x',x) 



.7=0 



C[J] = / V[9,9]e 



<e\A+j\9> 



(78) 
(79) 



A' 



Vc_ 

Vv 



E 



E-^k 2 -' 

2m 



l + v(v-Q)A> 



(83a) 



(83b) 



The fc-space summation can now be evaluated and the 
CPA equations solved. From these results the calculation 
of the density of states can be done. 
The effective medium operator that is defined by equa- 
tion (El) is 



A eS (k,k',z) 



2J.2 



2m 



Q\5ki 



(84) 



This result has already been obtained in an independent 
discussion of the CPA for electrons in a random potential 
[6(1 and can be interpreted as the continuum version of 
the classical lattice theories, e.g. Q. 
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